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Abstract.  The  paper  is  concerned  with  the  development  of  a  fast,  accurate,  and  versatile  method  of  simulat¬ 
ing  the  propagation  of  a  blast  wave  within  complex  built  environments.  An  ability  to  complete  a  simulation 
of  the  propagation  of  a  blast  wave  within  a  few  seconds  or  minutes  is  an  essential  tool  for  evaluating  its  im¬ 
pact  on  key  structures  and  to  find  an  optimal  design  for  components  such  as  blast  barriers.  Existing  methods 
of  modeling  the  propagation  of  a  blast  wave  fail  to  satisfy  all  of  the  above  requirements.  An  alternative  me¬ 
thod  proposed  here  is  to  use  a  simulation  approach  implemented  within  a  coarse  spatial  and  time  framework, 
where  the  mesh  elements  and  time  steps  are  orders  of  magnitude  larger  than  those  used  in  conventional  CFD 
simulations.  Preliminary  research  suggests  that  the  reduced  computational  load  resulting  from  this  method 
will  allow  simulations  to  be  executed  several  orders  of  magnitude  faster  than  conventional  CFD  methods, 
such  that  a  simulation  which  would  take  a  day  to  process  could  be  completed  within  a  few  seconds.  The  pa¬ 
per  describes  the  new  approach  in  detail. 
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Fluid  Dynamics,  very  fast  simulation. 


1  Introduction 

This  concept  paper  proposes  a  new  approach  to  simulating  the  propagation  of  blast  waves  in  a 
built-up  environment  that  has  the  potential  to  be  accurate,  generate  results  rapidly  (in  a  matter 
of  seconds),  and  has  the  versatility  to  model  complex  building  configurations.  Such  a  tool  is 
necessary  to  allow  engineers  to  design  new  structures  (including  new  buildings,  retrofits  of 
existing  bui ldings,  and  protective  structures  such  a s  b last  w alls)  that  pe rform  e ffectively  in 
terms  of  both  blast-mitigation  and  cost-effectiveness. 

Existing  blast  modeling  tools  are  a  trade-off  between  the  complexity  of  the  environment  they 
can  model  and  the  time  they  take  to  generate  results.  Modeling  tools  that  can  produce  results 
rapidly  are  direct  mapping  devices  (see,  for  example,  Remennikov  (2003)),  the  best  perform¬ 
ing  of  which  are  the  neural  network  models  such  as  those  described  by  Remennikov  and  Rose 
(2007)  and  Flood  et  al.  (2009).  Artificial  neural  networks  (an  introduction  to  which  is  pro¬ 
vided  by  Flood  &  Kartam  (1994))  are  very  versatile,  and  capable  of  considering  many  inde¬ 
pendent  variables  that  have  non-linear  relationships  with  the  dependent  variables.  The  Re¬ 
mennikov  a  nd  R  ose  ( 2007)  models  w  ere  t  rained  using  da  ta  from  m  iniature  bom  b-barrier- 
building  e  xperiments  ( Chapman  et  al.,  1995) ,  w  hile  t  he  F  lood  et  al.  ( 2009)  m  odels  w  ere 
trained  u  sing  da  ta  s  ynthesized  from  C  FD  ( Computational  F  luid  D  ynamic)  s  imulations  and 
other  established  modeling  techniques.  These  neural  network  models  can  produce  results  in  a 
fraction  of  a  second,  and  can  be  very  accurate.  F  igure  1 ,  for  example,  demonstrates  the  per¬ 
formance  of  a  neural  network  trained  to  predict  the  peak  pressure  (psi)  on  the  face  of  a  build¬ 
ing  where  there  is  a  blast  barrier  positioned  between  the  bomb  and  the  building,  for  the  set  of 
configurations  illustrated  in  Figure  2.  Each  point  in  the  figure  represents  1  of  252  randomly 
selected  test  problems,  plotting  the  neural  networks  predicted  value  against  the  actual  value  - 
if  the  model  was  perfect,  all  the  points  would  fall  on  the  diagonal  line.  The  results  gave  a  cor- 
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relation  between  the  model  output  and  actual  values  of  0.9959,  indicating  excellent  perfor¬ 
mance. 

Unfortunately,  empirically  derived  models  (such  as  neural  networks)  cannot  usually  extrapo¬ 
late  to  problems  beyond  those  represented  by  the  data  set  used  to  develop  the  model.  Moreo¬ 
ver,  the  size  of  the  data  set  required  to  develop  direct  mapping  empirical  models  (such  as  the 
above  artificial  neural  networks)  increases  geometrically  with  the  number  of  independent  va¬ 
riables  describing  the  problem.  I  n  practical  terms,  for  blast  wave  modeling  this  limits  the 
complexity  of  a  problem  to  about  five  independent  variables  -  this  has  constrained  applica¬ 
tion  to  setups  such  as  that  shown  in  Figure  2  comprising  a  two-dimensional  blast  wave  propa¬ 
gating  over  a  blast  barrier  onto  the  face  of  a  building,  where  the  barrier  and  building  are  per¬ 
pendicular  to  the  plane  of  the  blast  wave. 


252  Test  Patterns  Performance 


Fig.  1:  T ypical  P erformance  Results  for  an  Artificial  N eural  N etwork  Model  D eveloped  to 
Predict  Peak  Pressure  on  the  Surface  of  a  Building 
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Fig.  2:  Simple  Bomb-Barrier-Building  Configuration,  with  Blast  Wave  Acting  along  a  Ver¬ 
tical  Plane  Perpendicular  to  the  Blast  Barrier  and  Target  Building 

Blast  waves  propagating  through  more  complex  environments,  and  acting  in  three  spatial  di¬ 
mensions,  can  usually  only  be  considered  using  CFD  techniques  (such  as  AUTODYN 
(ANSYS,  2008)).  Unfortunately,  three  dimensional  CFD  models  of  blast  wave  propagation, 
even  when  limited  to  a  single  barrier  and  building  configuration  and  run  on  a  supercomputer, 
can  take  several  days  or  more  to  complete  a  single  simulation  run  (Flood  et  al,  2009). 

One  approach  for  simplifying  the  difficulty  in  blast  load  predictions  is  to  use  ray  tracing.  This 
uses  an  algorithm  that  identifies  the  most  significant  paths  (the  shortest)  that  a  blast  wave  can 
follow  from  the  point  of  detonation  to  specific  target  points,  taking  into  account  reflection  and 
diffraction.  The  time-based  forms  of  the  waves  arriving  along  each  path  are  determined  using 
the  semi-empirical  TNT  Standard  Methodology  and  are  then  superimposed  using  the  LAMB 
Shock  Addition  rules  (Needham  &  Crepeau,  1981).  Enhancements  to  the  approach  have  been 
used  by  Frank  et  al.  (2007a,  2007b)  to  predict  behavior  in  environments  with  complex  geome¬ 
tries.  The  approach  certainly  provides  a  highly  versatile  method  of  modeling  complex  inter¬ 
nal  geometries,  and  it  is  claimed  that  the  results  are  of  reasonable  accuracy  and  that  the  model 
runs  fast.  However,  a  comprehensive  set  of  case  studies  is  required  to  determine  more  pre¬ 
cisely  the  accuracy  and  processing  speed  relative  to  established  approaches  to  the  problem. 
The  algorithm  required  to  determine  all  significant  paths  for  the  blast  wave  appears  to  be  too 
complex  to  allow  results  to  be  generated  in  a  matter  of  seconds  for  all  target  points  across  all 
relevant  surfaces  of  the  environment  that  would  be  required  for  the  applications  proposed  in 
this  paper. 

An  alternative  approach  to  these  issues,  considered  by  Lohner  et  al.  (2004),  was  to  test  the 
sensitivity  of  processing  time  and  accuracy  on  the  coarseness  of  the  modeling  mesh  for  three 
dimensional  CFD  simulations  of  blast  wave  propagation  through  complex  building  geome- 
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tries.  In  an  example  study  of  a  concert  hall,  consideration  was  given  to  a  range  of  resolutions 
ranging  from  main  element  sizes  of  0.3  m  to  1 .2  m  in  length.  It  was  found  that  moving  to  the 
coarser  mesh  reduced  processing  time  from  18  hours  to  7  minutes,  although  the  predictions  of 
the  coarse  mesh  model  were  about  50%  off  compared  to  the  fine  mesh  model.  While  the 
speed  o  f  p  rocessing  o  f  t  he  co  arse  m  esh  ap  proach  m  akes  it  a  ccessible  t  o  u  sers  o  f  d  esk-top 
computers,  the  authors  of  this  paper  believe  greater  accuracy  in  the  predictions  of  the  model  is 
needed. 


2  Coarse-Grain  Based  Simulation 

This  paper  proposes  a  coarse-grain  approach  to  achieving  a  modeling  system  that  is  fast,  accu¬ 
rate,  and  versatile.  It  differs  fundamentally  from  the  coarse-grain  approach  of  Ldhner  et  al. 
(2004)  discussed  above  in  that  it  uses  empirical  rather  than  theoretically  derived  functions  to 
drive  a  simulation.  In  Lohner’s  study,  the  coarseness  of  a  model  was  achieved  by  simply  in¬ 
creasing  the  size  of  the  spatial  elements  comprising  a  model,  while  using  the  same  discretized 
driving  equations  used  in  the  fine-grain  models.  Increasing  the  size  of  the  elements  reduces 
dramatically  the  number  required  for  any  given  situation  and  thus  similarly  reduces  the  com¬ 
putational  load  of  a  simulation,  hence  the  significant  reduction  in  processing  time.  However, 
the  driving  equations  used  assume  an  infinitesimally  small  element  size,  and  so  become  less 
accurate  the  larger  the  size  of  those  elements  -  this  is  also  true  of  the  size  of  the  time  step  dur¬ 
ing  the  execution  of  a  simulation.  Consequently,  the  accuracy  of  the  simulations  reduces  very 
quickly  with  model  coarseness.  In  this  study,  it  is  proposed  to  compensate  for  this  in  the  fol¬ 
lowing  two  ways: 

1.  Artificial  neural  networks  will  be  used  to  leam  a  more  accurate  form  for  the  driving 
equations  operating  in  a  coarse-grain  modeling  environment.  These  neural  networks 
will  be  trained  based  on  data  gathered  from  a  comprehensive  set  of  CFD  simulations 
of  building  element  geometries. 

2.  Each  co  arse-grain  e  lement  w  ill  in  elude  a  n  additional  s  et  o  f  independent  v  ariables 
sampling  the  state  of  the  element  historically  (in  the  time  domain)  to  compensate  for 
the  loss  of  information  by  the  low  resolution  in  the  spatial  domain. 

This  a  pproach,  d  eveloping  a  s  et  o  f  d  iscretized  driving  e  quations  tailored  t  o  a  co  arse-grain 
modeling  environment,  has  been  demonstrated  to  be  very  effective  in  an  earlier  series  of  stu¬ 
dies  modeling  dynamic  heat  transfer  in  complex  building  configurations  (Flood  et  al.  (2004)). 
For  transient  heat  transfer,  it  was  found  that  a  system  could  be  decomposed  into  very  coarse 
elements  w  hile  m  aintaining  accu  racy  in  pr  edicted  pe  rformance.  I  ndeed,  e  ach  w  all,  c  eiling, 
floor,  or  open  space  could  be  represented  by  a  single  element,  as  shown  in  F igure  3.  T  wo 
primary  element  types  were  considered,  those  representing  open  spaces  (such  as  rooms  and 
attic  spaces)  and  those  representing  boundaries  between  spaces  (such  as  walls,  floors/ceilings, 
and  roof  membranes).  The  simulation  advances  in  a  manner  very  similar  to  CFD  techniques, 
by  predicting  the  temperature  at  each  element  at  the  next  point  in  time  based  on  its  current  and 
historic  t  emperatures  a  nd  th  ose  o  f  its  n  eighboring  e  lements  -  the  d  ifference  b  eing  th  at  th  e 
equations  used  to  make  these  predictions  are  empirically  derived  and  compensate  for  the  large 
distances  between  adjacent  points  and  between  points  in  time  by  using  many  historical  sam¬ 
pling  points.  In  a  case  study  of  a  simple  house,  simulated  for  a  one-year  heating  and  cooling 
cycle  with  a  1 5  minute  time-step,  the  model  was  found  to  end  within  0.25  °F  of  the  actual  val¬ 
ue.  T  he  results  also  indicated  that  a  three-dimensional  coarse-grain  simulation  of  the  house 
could  be  performed  several  thousand  times  faster  than  a  conventional  CFD  model. 
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(g)  =  boundary  modeling  element 
O  =  space  modeling  element 


o 


0  O  0  O  0 

— 0 — I - 0 - 


Fig.  3:  Example  Decomposition  of  a  Simple  House  into  Coarse-Grain  Elements  for  Modeling 
Dynamic  Heat  Transfer 


3  Coarse-Grain  Modeling  Applied  to  Blast  Wave  Propagation 

There  is  a  critical  difference  between  modeling  heat  transfer  and  blast  wave  propagation  that 
invalidates  the  coarse  grain  approach  described  in  S  ection  2 .  Fundamentally,  the  problem 
arises  from  the  fact  that  the  temperature  distribution  within  a  structure  changes  very  gradually 
in  the  spatial  domain  relative  to  the  distance  between  the  modeling  elements,  whereas  a  blast 
wave  has  a  v  ery  short  wavelength  compared  to  the  distance  between  coarse  grain  elements. 
Consequently,  as  illustrated  in  Figure  4,  an  advancing  blast  wave  may  be  lost  within  the  mesh, 
perhaps  never  intercepting  more  than  one  coarse  grain  element  at  a  time,  and  for  much  of  its 
existence  not  intercepting  any  coarse  grain  elements.  T  his  means  that  there  would  never  be 
enough  information  about  the  state  of  the  blast  wave  to  be  able  to  make  predictions  about  its 
state  at  a  succeeding  point  in  time.  The  coarse  grain  approach  as  described  for  modeling  tran¬ 
sient  heat  transfer  is  clearly  inadequate. 
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Fig.  4  :  Illustration  of  the  Fact  that  an  Advancing  Blast  Wave  will  Intercept  Very  Few  Coarse 
Grain  Elements  at  any  Point  in  Time 

The  basis  of  the  proposed  solution  to  this  problem  is  to  identify  the  state  of  the  blast  wave  by 
all  co arse-grain  elements  adjacent  to  it,  not  just  those  intercepted  by  it,  and  to  advance  the 
simulation  at  each  step  by  jumping  to  the  time  at  which  the  blast  wave  intercepts  the  next 
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coarse  grain  element.  The  key  properties  of  the  wave  to  be  captured  at  any  point  in  time  dur¬ 
ing  t  his  pr  ocess  are  i  ts  1  ocation,  c  orresponding  peak  pressure,  t  rough  pressure,  a  nd  wave 
length,  as  indicated  in  Figure  5. 


Fig  5:  Properties  of  Blast  Wave  to  be  Captured  by  the  Proposed  Coarse  Grain  Modeling  Ap¬ 
proach 

The  s  pecific  s  teps  i  n  t  he  pr  oposed  s  imulation  c  an  b  e  unde  rstood  by  r  eferring  t  o  F  igure  6 . 
Here,  the  advancing  blast  wave  is  represented  by  the  gray  line,  moving  south  west,  and  is  reg¬ 
istered  at  all  adjacent  coarse  grain  elements  shown  as  color-shaded  active  nodes.  T  he  blue 
nodes  are  trailing  the  blast  wave,  the  yellow  one  is  currently  intercepted  by  the  wave,  and  the 
red  nodes  are  leading  the  wave.  Each  of  these  nodes  registers  the  following:  the  time  until  (or 
since)  the  peak  pressure  of  the  blast  wave  reaches  (or  reached)  that  node;  the  peak  and  trough 
pressures  that  will  occur  (or  occurred)  at  that  node;  and  the  velocity  of  the  propagating  wave 
(presented  in  rectangular  vector  components  so  that  information  on  the  direction  of  travel  of 
the  wave  is  also  implied).  To  advance  the  simulation  to  the  next  point  in  time,  the  following 
steps  (which  similar  to  next-event  based  discrete-event  simulation)  must  be  implemented: 

STEP  1 :  Find  the  leading  active  node  (red  in  this  figure)  that  has  the  shortest  time  until  ar¬ 
rival  of  the  blast  wave. 

STEP  2 :  Advance  the  simulation  to  this  point  in  time  by  increasing  by  a  corresponding 
amount  the  time  since  arrival  at  the  trailing  active  nodes,  and  reducing  the  time  before  ar¬ 
rival  of  the  leading  active  nodes. 

STEP  3:  Remove  any  nodes  from  the  active  list  that  are  no  1  onger  adjacent  to  the  blast 
wave. 

STEP  4 :  Add  all  nodes  to  the  active  list  that  are  now  adjacent  to  the  blast  wave.  T hese 
nodes  will  be  leading  the  blast  wave. 
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STEP  5:  For  each  leading  node  that  has  just  been  added  to  the  active  list  in  STEP  4,  pre¬ 
dict  the  time  until  arrival  of  the  blast  wave,  the  peak  and  trough  pressures  on  its  arrival, 
and  the  velocity  the  blast  wave  will  have  on  arrival  at  that  node  (in  rectangular  vector 
coordinates).  These  predictions  must  be  based  on  the  state  of  the  blast  wave  registered  at 
its  neighboring  nodes.  In  Figure  6  for  example,  the  green  arrows  point  to  a  leading  node 
that  has  just  become  active,  and  emanate  from  the  nodes  that  will  be  used  to  calculate  its 
blast  wave  attributes.  In  the  proposed  approach,  these  predictions  will  be  made  by  artifi¬ 
cial  neural  network  modules  trained  using  data  from  specific  CFD  experiments. 


Key: 

o  Trailingactive  nodes 
o  Current  active  nodes 
•  Leading  active  nodes 


J 


coarse  grain 


elements 


Fig  6:  An  Advancing  Blast  Wave  and  its  Adjacent  Active  Nodes  (Coarse  Grain  Elements) 


4  Estimate  of  Performance  of  Proposed  Coarse  Grain  Modeling  Approach 

The  next  stage  in  the  project  will  be  concerned  with  implementing  the  proposed  method  and 
evaluating  its  performance  in  terms  of  processing  speed  and  accuracy.  However,  based  on 
past  experiences  w  ith  n  eural  ne  twork  m  ode  ling  a  nd  c  oarse  g  rain  m  odeling  i  t  i  s  pos  sible  t  o 
make  a  rough  estimate  of  likely  processing  speed  performance.  This  will  be  completed  with 
reference  to  the  2  di  mensional  coarse  grain  model  shown  in  F igure  7,  that  includes  a  7  m 
square  grid  with  coarse  grain  elements  located  at  1  m  spacings,  an  obstruction  on  the  top-left 
side,  and  a  charge  placed  at  the  hot  tom- left  comer.  The  performance  of  this  model  will  be 
compared  t  o  t  hat  o  f  a  conventional  C  FD  m  odel  ( that  us  es  the  F  inite  D  ifference  M  ethod  - 
FDM)  that  places  elements  on  a  2  cm  mesh  (the  largest  spacing  found  in  previous  experi¬ 
ments  to  provide  an  acceptable  level  of  accuracy). 
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Fig  7:  Example  Two-Dimensional  Coarse  Grain  Model 

A  summary  of  the  processing  load  for  both  the  coarse  grain  and  CFD  models  of  this  scenario 
are  shown  in  Table  1.  At  each  step  in  the  coarse  grain  simulation  only  one  element  on  aver¬ 
age  will  be  recalculated,  whereas  the  CFD  model  will  have  approximately  122,500  elements 
to  recalculate  (see  column  2  i  n  T  able  1) .  In  addition,  whereas  the  coarse  grain  model  will 
require  71  time  steps  to  advance  the  blast  wave  from  the  bomb  to  the  top  right  element  (re¬ 
membering  there  will  be  additional  time  steps  resulting  from  the  fact  that  the  wave  will  be 
reflected  of  f  t  he  obs  traction),  experience  i  ndicates  t  hat  the  C  FD  s  imulation  w  ill  r  equire 
around  7,425  time  steps  to  make  the  same  advancement.  H  owever,  a  d  isadvantage  for  the 
coarse  grain  approach  is  that  the  processing  load  to  update  one  element  at  any  time  step  may 
be  around  600  times  greater  than  the  CFD  model.  This  is  because  the  CFD  model  uses  a  very 
simple  discretized  function  whereas  the  coarse  grain  model  requires  the  use  of  an  artificial 
neural  network  -  the  value  of  600  i  s  based  on  earlier  experiments  with  neural  networks  to 
model  the  state  transition  function  for  a  coarse  grain  simulation.  The  net  result  of  all  these 
factors  is  a  processing  speed  advantage  of  about  21,351  for  the  coarse  grain  approach.  In  oth¬ 
er  words,  if  the  CFD  simulation  took  1  da  y  to  process,  the  coarse  grain  model  would  take 
about  4  seconds.  This  advantage  would  be  expected  to  increase  for  a  three  dimensional  mod¬ 
el.  For  example,  using  the  same  2  cm  mesh  size,  the  CFD  model  would  require  50  layers  of 
elements  for  every  meter  depth  whereas  the  coarse  grain  model  would  require  just  one  layer. 
This  would  increase  the  processing  advantage  by  a  factor  of  50  so  that  a  CFD  simulation  that 
took  1  year  to  process  could  be  completed  by  the  coarse  grain  model  in  just  30  seconds. 


Table  1:  Processing  Speed  Comparison  of  Conventional  CFD  Model  (Finite  Difference 
Method  -  FDM)  and  Proposed  Coarse  Grain  Method 


Model 

Type 

Number 
Elements/ 
Time  Step 

Time 
Steps 
(to  far 
corner) 

Processing 
Load  Per 
Calculation 
(normalized) 

CFD 

Equivalent 

Calculations 

Processing 

Speed 

Ratio 

CFD 

1  Day 
Processing 

Extrapolate 
to  3D  model 
Processing 
Speed  Ratio 

Extrapolate 
to  3D  model 

1  year 
processing 

CFD 

122,500 

(2cm  mesh) 

7,425 

1 

909,562,500 

1 

1  day 

1 

1  year 

Coarse 

Grain 

1 

71 

600 

(6x100) 

42,600 

21,351 

4  seconds 

1,067,550 
(2D  x  50) 

30  seconds 

8 


Such  an  increase  in  performance  would  certainly  make  three  dimensional  simulation  of  blast 
wave  propagation  in  a  complex  built  environment  accessible.  However,  the  other  question  to 
be  answered  is  the  accuracy  of  the  predictions.  This  is  a  more  difficult  question  to  address 
than  that  of  processing  speed  without  actually  performing  a  comprehensive  set  of  customized 
experiments.  Such  a  set  of  experiments  is  proposed  for  the  next  stage  of  this  study 


5  Discussion 

Empirical  methods  of  modeling  the  effects  of  bomb  blast  waves  on  buildings  are  fast  and  of¬ 
ten  accurate,  but  lack  the  flexibility  to  consider  anything  other  than  the  simplest  of  built  envi¬ 
ronment  geometries.  On  the  other  hand,  the  CFD  models,  while  potentially  very  versatile,  can 
take  several  days  to  process  a  three  dimensional  model.  As  a  consequence,  CFD  models  can¬ 
not  be  used  for  interactive  design,  do  not  allow  engineers  to  consider  more  than  a  few  alterna¬ 
tive  setups,  and  prohibit  consideration  of  stochastic  effects  (such  as  uncertainty  about  the  size 
and  location  of  the  bomb,  the  open/closed  status  of  blast  doors,  etc.)  since  this  would  require 
multiple  simulation  executions  using  techniques  such  as  Monte  Carlo  sampling. 

In  response,  this  paper  has  proposed  a  new  method  of  modeling  based  on  simulating  with  a 
coarse-grained  modeling  mesh.  Normally  coarse  meshes  result  in  significant  reductions  in  the 
accuracy  of  predictions;  the  intent  here,  however,  is  to  compensate  for  the  loss  of  spatial  in¬ 
formation  by  developing  custom  functions  for  coarse  grain  analysis  that  also  sample  in  the 
time  domain.  P  revious  studies  in  the  field  of  transient  heat  transfer  have  shown  this  to  be 
viable  if  artificial  neural  networks  are  used  to  learn  the  functions  that  drive  the  simulation. 
The  approach  has  the  potential  to  increase  processing  speed  by  several  orders  of  magnitude 
compared  to  conventional  CFD  simulations.  The  next  step  in  this  project  is  to  implement  the 
new  method  and  validate  it  for  a  comprehensive  range  of  complex  building  environments. 
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